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We present a detailed investigation of the magnetic and structural properties of magnetically 
doped 3D topological insulator Bi2Se3. From muon spin relaxation measurements in zero magnetic 
field, we find that even 5% Fe doping on the Bi site turns the full volume of the sample magnetic 
at temperatures as high as ~ 250K. This is also confirmed by magnetization measurements. Two 
magnetic "phases" are identified; the first is observed between ~ 10 — 250K while the second appears 
below ~ 10K. These cannot be attributed to impurity phases in the samples. We discuss the nature 
and details of the observed magnetism and its dependence on doping level. 
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I. INTRODUCTION 

Topological insulators (TIs) are 3D materials with 
a band gap in the bulk-electronic structure but with 
gapless, delocalized, surface states^. The surface-states 
are protected by time-reversal symmetry (TRS), and 
cannot be gapped without changing substantially the 
bulk-electronic structure. This special relation between 
the bulk and surface states is the basic difference be- 
tween these systems and systems showing ordinary sur- 
face states. The protected surface states are believed 
to be robust to disorder, interactions, and thermal fluc- 
tuations, potentially leading to room-temperature de- 
vice applications^. In recent years, doped TIs have 
been investigated due to potential new phenomena when 
the topological surface states interact with impurities or 
other electronic states in the bulk. For example, it was 
found that by intercalating a prototypical TI, Bi2Se3, 
with Cu it becomes superconducting below ~ 4KP. More 
recently, the interaction between the surface states of a 
TI and magnetic impurities have also been attracting at- 
tention since such impurities are expected to break TRS. 
In this context, an important issue is the stability of the 
TI surface states to the presence of impurities. 

In the bulk, it was shown that a TI with magnetic 
impurit ies can have long range magnetic order in the 
metallicPE' and insulating^ phases. On the surface, long 
range order can also be formed independent of the bulk 
magnetic ordering, as the RKKY interaction induced by 
the Dirac fermions is generally ferromagnetic when the 
Fermi energy (Ep) is close to the Dirac point. Both ef- 
fects can lead to breaking of TRS, resulting in a gap 
opening at the Dirac point that makes the surface Dirac 
fermions massive. Chen et alPhave shown that the Dirac 
gap can be observed in magnetically doped samples with 
or without bulk ferromagnetism. Furthermore, it was 
found that Ep can be tuned into the surface-state gap 
which in turn may allow the observation of many interest- 



ing topological phenomena, such as the image magnetic 
monopole induced by an electric point charge 89 . The- 
oretical work on magnetically doped TIs has also been 
very active^HSl. These studies predict a variety of effects 
due to the RKKY interactions between different mag- 
netic impurities mediated by the surface states. For ex- 
ample, Zhu et al. 1 - predict that such interaction consists 
of Heisenberg-like, Ising-like, and Dzyaloshinskii-Moriya- 
like terms. They devise a new way to control surface 
magnetism using electrical field or even by changing the 
average spacing between different impurities at the sur- 
face, e.g., by changing their doping level. However, most 
importantly this surface magnetism is predicted to be 
dramatically different from the bulk due to the different 
type of magnetic interactions at the surface. 



In this paper we investigate magnetically (Fe) doped 
Bi 2 Se 3 , at different doping levels. Using zero field (ZF) 
muon spin relaxation (/iSR) and bulk magnetization mea- 
surements, we find that even at low doping levels, the 
full volume of the material becomes magnetic at a rela- 
tively high temperature (~ 250K). We identify two mag- 
netic regimes; the first between ~ 10 — 250K and the 
second below ~ 10K. The nature of the magnetism in 
both regimes, as well as the effective size of the magnetic 
moment per doping Fe, are independent of the doping 
level. However, the dynamics due to magnetic fluctua- 
tions increase at lower doping. Moreover, near ~ 250K, 
the size of local magnetic field at the muon site increases 
more sharply as the temperature is decreased for higher 
Fe doping. These results indicate that Fe-Fe interactions 
dominate the magnetism at high temperatures, while the 
single Fe ion properties dominate in the low temperature 
regime. These studies form a basis for future studies of 
the magnet ism n ear the surface using depth resolved low 
energy /^SR [ 15 * 16 \ which may shed light on the interplay 
between magnetism and topological surfaces states in this 
important class of materials. 



II. EXPERIMENTAL 

Single crystals of Bi 2 _a;Fe x Se3, with x = 0, 0.1, 0.2 
and 0.3, were grown from a melt using the Bridgman 
method. Corresponding amounts of Bi, Fe, Se of mini- 
mum purity 99.99% were mixed and sealed in evacuated 
silica ampules. The ampules were annealed at 820° C 
over 10 h for homogenization. The melt was then cooled 
down to 640° C at a rate of 4° C/h and then quenched 
into cold water. Well-formed silvery crystal rods were 
obtained, which could be easily cleaved into plates with 
flat shiny surfaces. The magnetization measurement re- 
ported here were obtained using a Physical Properties 
Measurement System (PPMS, Quantum Design). 

The Bi2_ a; Fe a ;Se3 (x = 0.2, 0.3) crystals were charac- 
terized at room temperature by powder X-ray diffraction 
(XRD) using a D8 Advance Bruker AXS diffractometer 
with Cu K Q radiation. Additionally, the Bi 2 _ a; Fe a ;Se3 
(x = 0.2) crystal was studied by means of neutron pow- 
der diffraction (NPD) at the SINQ spallation source 
at the Paul Scherrer Institute (PSI, Switzerland), us- 
ing the high-resolution diffractometer for thermal neu- 
trons HRPT 1 ? (wavelengths A = 1.494 A and 1.155 A). 
The NPD patterns were measured at temperatures in the 
range 2 — 300 K. For the diffraction measurements, part 
of the crystal was powderized and loaded into a vana- 
dium container with an indium seal in a He glove box. 
Refinement of the crystal structure parameters was done 
using FULLPR0F software^, with its internal tables for 
neutron scattering lengths. 

The fj.SK experiments were performed on the DOLLY 
spectrometer at PSI, Switzerland. In these experiments 
100% polarized (along the beam direction, z) positive 
muons are implanted in the sample. In our measure- 
ments single crystals were suspended on an aluminizcd 
Mylar tape and mounted into a He gas flow cryostat. The 
crystals were aligned such that the initial polarization of 
implanted muons is along the c-axis (and the applied field 
in non-zero field measurements). Each implanted muon 
decays (lifetime 7), = 2.2 /isec) emitting a positron pref- 
erentially in the direction of its polarization at the time 
of decay. Using appropriately positioned detectors, one 
measures the asymmetry of the muon beta decay along 
z as a function of time A(t), which is proportional to the 
time evolution of the muon spin polarization. A(t) de- 
pends on the distribution of internal magnetic fields and 
their temporal fluctuations. 



III. RESULTS 

A. Magnetization measurements 

In Fig. [T] we show the magnetic susceptibility, x, as a 
function of temperature. In this figure we present both 
zero field cooled (ZFC) and field cooled (FC) measure- 
ments in a magnetic field of B = 100 mT applied per- 
pendicular to the crystallographic c-axis (B 1 c). The 
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FIG. 1: The FC (blue) and ZFC (red) susceptibility mea- 
sured in B = 100 mT magnetic field (B _L c) as a function of 
temperature for the different Fe doping. The solid lines are 
a guide to the eye and the insets are a zoom of the shaded 
areas (0-30 K). 



mother compound Bi2Se3 gives a small diamagnetic sig- 
nal with almost no temperature dependence. In contrast, 
Fe doping makes Bi2- x Fe 2; Se3 magnetic below ~ 250K, 
confirmed by a clear split between the ZFC and FC sus- 
ceptibility. The magnetic nature seems to change at 
temperatures below ~ 10 K, where a small abrupt up- 
turn is observed for all doping levels measured (see in- 
sets of Fig. [I]). Similar behaviour is observed in mea- 
surements with the magnetic field applied parallel to the 
c-axis (B |j c), but with a more pronounced upturn below 
~ 10 and smaller x values in general. 

Measurements of the hysteresis loops in all four sam- 
ples are shown in Figs. | (B 1 c) and Fig. [3] (B | c). 
These were measured at two different temperatures; 2 
and 35 K to better show the difference above and below 
10 K. The diamagnetic and temperature independent na- 
ture of the mother compound Bi2Se3 (x = 0) is evident 
in the absence of hysteresis and the negative slope of M 
vs. B yielding x ~ — l -6 emu/g. However, for all doped 
samples we observe a clear opening of the hysteresis when 
the field is applied perpendicular to the c-axis (Fig. [2]) . In 
contrast, we find that the hysteresis loops are much nar- 
rower when the field is applied along the c-axis (Fig. |3| . 
This strong anisotropy indicates that the magnetic easy 
axis in Bi2- I Fe x Se3 lies within the Bi layer and almost 
perpendicular to the c-axis. Furthermore, the hysteresis 
loops seem to be narrower below 10 K pointing to the 
possibility of a different orientation of the magnetic easy 
axis between the high and low temperature regimes. 

Note, in the Blc measurements (Fig. [2| the rema- 
nence depends only weakly on temperature, but the satu- 
ration of the magnetization differs significantly. In addi- 
tion, we find that the magnetization saturates gradually 
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FIG. 2: The magnetization as a function of magnetic field 
(B _L c) for the different Fe doping. The red squares and blue 
circles are measurements at T — 2 and 35 K, respectively. 
The solid lines are a guide to the eye. 




FIG. 3: The magnetization as a function of magnetic field 
(B || c) for the different Fe doping. The red squares and blue 
circles are measurements at T — 2 and 35 K, respectively. 
The solid lines are a guide to the eye. 



at 2K, while at 35K it is much clearer, though lower in 
magnitude. This difference is again a possible indica- 
tion of the different nature of magnetic ordering below 
10 K, e.g. due to a slight reorientation of the magnetic 
easy axis. Our measurements are similar to those re- 
ported earlier^ on similar samples, though the magne- 
tization values measured in our samples are significantly 
higher for similar doping levels. As we show below, XRD, 
NPD and /iSR measurements strongly indicate that the 
observed magnetization cannot be attributed to ferro- 



magnetic inclusion of Fe-Se binary compounds. 



B. NPD and XRD measurements 

Both room-temperature NPD and XRD experiments 
show that the samples consist mainly of the rombohe- 
dral Bi2_ 2; Fe a: Se3 phase (space group R-3m). However, 
we also detect a small impurity phase of Fe 3 Se4 (about 
5% mass), which has a monoclinic structure and orders 
ferrirnagnetically below T c ~ 314 K 19 . Difference Fourier 
density map analysis reveals that Fe substitutes Bi atoms 
in the structure. The refined Fe occupancy is in agree- 
ment with the nominal composition, but the deduced val- 
ues contain large systematic errors. This is due to the 
presence of the small impurity phase and a strong pre- 
ferred orientation of the flat-shape crystallites in the pow- 
der along (OOl)-direction. This in turn affects their pack- 
ing and the resulting diffraction patterns seen in both 
XRD and NPD. We account for this artifact by using 
the March formula in the FULLPR0F refinement analysis. 
Note, since the scattering lengths of Bi and Fe are quite 
close (8.53 fm and 9.45 fm, respectively), NPD is not 
very sensitive to the Bi/Fe-occupancy ratio. Therefore, 
their sum was simply fixed at 2. High statistics NPD 
data, which were collected at two different temperatures 
2K and 20K, do not exhibit any statistically significant 
difference due to possible magnetic contribution. This 
may be due to a small magnetic moment on the Fe ions 
or due to a disordered nature of the magnetic phase in 
these compounds. The crystal structure parameters for 
the main phase refined from XRD (300K) and NPD (2K 
and 300K) are presented in the Table [I] 



C. /iSR measurements 

Compared to conventional dc magnetization measure- 
ments, /iSR is a much more sensitive method to inves- 
tigate the magnetism of these compounds. It also has 
two very important advantages: (I) detecting intrinsic 
magnetism even in zero applied magnetic fields and (II) 
since it is a local probe method, the contribution of im- 
purity phases is proportional to their volume fraction and 
can be easily disentangled from the contribution of the 
studied compound. Typical /iSR asymmetries measured 
in the x = 0.2 sample are shown in Fig. [4] Similar be- 
haviour was observed in the other two Fe doped com- 
pounds. Note, in zero applied field, the asymmetry re- 
laxes to finite value. A priori, one cannot determine the 
source of this relaxation. It can be due to either a dis- 
tribution of internal static fields, dynamic fields, or a 
combination of the two. However, it is evident in Fig. HI 
that application of an external magnetic field along the 
initial direction of the muon spin (so-called longitudinal 
field or LF-^SR) decouples the relaxation^. This decou- 
pling, together with the appearance of "wiggles" in the 
relaxation curves when the field is applied, are evidence 





XRD (300K) 


NPD (300K) 


NPD (2K) 


a (A) 


4.14072(10) 


4.13914(4) 


4.12518(3) 


c(A) 


28.64692(60) 


28.63151(56) 


28.47718(40) 


Bi (0,0,2) 


0.40057(6) 


0.40107(14) 


0.40125(9) 


B Bi (A 2 ) 


0.997(52) 


1.743(67) 


0.498(45) 


Se2 (0,0,2) 


0.20984(10) 


0.21071(8) 


0.21062(6) 


B So2 (A 2 ) 


3.173(144) 


1.393(97) 


0.745(65) 


l~Lp i-Tl"wp j-ttexp 


6.55, 8.46, 3.77 


3.80, 4.75, 3.65 


3.46, 4.46, 2.00 


x 2 


5.04 


1.69 


4.99 



TABLE I: Crystal structure parameters of Bi 2 -^Fe ;c Se3 (x = 0.2) obtained from XRD (300 K) and NPD (300 and 2 K) refined 
in the space group R-3m (no. 166) with Sel in (3a) position (0, 0, 0), Bi and Fe in (6c) position (0, 0, z), and Se2 in (6c) position 
(0, 0, z). Bsi and Bs e 2 are the atomic displacement parameters of Bi and Se2 atoms, respectrively. R p ,R wp and R exp are the 
reliability factors and x 2 is the global chi-square (Bragg contribution). 
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FIG. 4: The asymmetry measured in the x = 0.2 compound 
with different applied longitudinal magnetic fields and tem- 
peratures: (a) 35K, (b) 5K and (c) 1.6K. The lines are fits to 
Eq. 0. 

of the quasi-static nature of the local fields experienced 
by the implanted muons, i.e., the local field is predomi- 
nantly static with a smaller dynamic component. In this 
case the relaxation can be described using a phenomeno- 
logical function; a Lorentzian Kubo-Toyabe (LKT) mul- 
tiplied by an exponential 21 . The LKT accounts for the 
relaxation due to the distribution of static fields while the 
exponential accounts for the depolarization due to the 
small dynamic component. The asymmetries in Fig. HI 
were fit to a function of the form 



A(i) = A P LK T(S,A,i)e-* /Tl , 



(1) 



where Aq is the initial asymmetry, Plkt{B, A, t) is the 
depolarization function due to a Lorentzian distribution 
of static fields with width A in the applied field B, and 
1/Ti is the spin lattice relaxation rate due to the fluctu- 
ating transverse magnetic field component at the muon 
site. In order to obtain a reliable fit in Fig. |1J we assume 
a temperature and field independent value of Aq and a 
field independent (but temperature dependent) value of 
A, with only 1/Ti allowed to vary as a function of the 



applied field. These assumptions produce reliable fits 
and very good agreement with the measured asymme- 
tries. Note that Eq. account for the full signal, and 
therefore no impurity phase (e.g. FesSe,^ contribution is 
detected in our measurements. This could be due to the 
small volume fraction of this phase and its correspond- 
ingly small contribution to our measured signal. Alter- 
natively, the single crystals used for these measurements, 
which were cleaved from large crystal rods, are free from 
impurities. 

Typically, the muon spin depolarization function due 
to a randomly oriented quasi-static local magnetic field 
exhibits a dip at early times, followed by a recovery to 
~ 1/3 of the initial asymmetry which continues to relax 
in an exponential-like manner at longer times. The ab- 
sence of a clear dip in the ZF measurements is due to the 
relatively small value of A (~ 1MHz) and the dynamic 
component of the local field. This makes it difficult to fit 
reliably all our data to Eq. (fTl). Instead we parametrize 
the temperature dependence in all three Fe doped com- 
pounds using a function of the form, 



A(t) = A 



,-xt 



1 



,-t/Ti 



(2) 



Here, A is a parameter proportional to A (approximately 
A ~ 2A) and quantifies the width of the static field dis- 
tribution. From these fits we extract the values of A and 
1/Ti as a function of temperature shown in Fig.plfor the 
three Fe doped samples. Note that, in agreement with 
the magnetization measurements, we observe a large in- 
crease in A below ~ 250 and ~ 10K. These coincide with 
broad and weak peaks in the values of 1/Ti, indicating 
that the Fe doped Bi2Se3 undergoes two different types 
of magnetic ordering at these temperatures. However, 
while A below ~ 250 K seems to depend strongly on the 
doping level, it is almost independent of x below ~ 10 K. 
In fact, the main difference between the doping levels ap- 
pears just below 250 K, where the increase in A is sharper 
for higher doping. The spin lattice relaxation rate, 1/Ti, 
depends more strongly on x. Its value increases with de- 
creasing x, indicating that the fluctuations of local fields 
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peratures shows that it increases linearly with x (Fig. [6J) . 
This behaviour indicates that the size of the magnetic 
moment of Fe does not depend strongly on x and that 
the increase in magnetization simply scales with the num- 
ber of Fe ions in the sample. In fact, from the slope of 
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FIG. 6: The saturation magnetization as a function of x, the 
Fe doping level at 35K and 2 K. The dashed lines are a linear 
fit of the data. 



FIG. 5: The temperature dependence of A (bottom) and 1/Ti 
(top) obtained from the ZF pSR data for the different Fe 
doping. 



are faster for lower doping. In all three doping levels 
a sharp peak is observed at ~ 5 K, which we attribute 
to the magnetic phase transition below ~ 10 K. In con- 
trast, the peak at high temperature is broader for all x 
values and it seems to shift to lower temperatures for 
larger Fe doping. We attribute this peak to the magnetic 
phase transition occurring below ~ 250 K. The width 
and quality of the results do not allow accurate determi- 
nation of the exact transition temperature, nevertheless, 
they indicate a more disordered and possibly inhomoge- 
neous magnetic transition (typically seen in dilute spin 
glasse:P2l). 

In contrast to the doped samples, the muon spin relax- 
ation rate in the mother compound is very small and has 
a very weak temperature dependence. Fitting the relax- 
ation to an exponential function gives a rate of ~ 0.008 
/is -1 (bottom panel of Fig. p3j). The only source of relax- 
ation in this case is expected to be dominated by the Bi 
nuclear magnetic moments. However, the observed relax- 
ation is much smaller than one expects from dipolar fields 
of nuclear moments and therefore strongly indicates that 
the muons stopping sites are far away from Bi. In our 
case, the muons are most probably located in the van der 
Waals gap between layers of Bi 2 Se 3 . 



IV. DISCUSSION 

We start by discussing the bulk magnetization results. 
The values of saturation magnetization, M sat , extracted 
from Fig.[2j are plotted as a function of x in Fig. [6] For all 
doping levels we find that M sat is higher at T = 2 K than 
35 K. Moreover, a close inspection of M sat at both tem- 



M sat as a function of x we find that the contribution of 
each Fe ion to the magnetization is equal (within error) 
at both temperatures and amounts to 1.91 ±0.06 emu/g. 
From M sat we estimate a magnetic moment per Fe of 
0.34(7) and 0.25(4) fi B at 2 and 35 K, respectively. Note, 
however, this contribution must decrease for much lower 
doping in order to approach the non-magnetic x — > 
limit. 

The fiSK results, summarized in Fig. [5] show that 
the width of the static field distribution (or A), which 
is proportional to the local magnetization, depends only 
weakly on x below ~ 50 K. Since the muon is a local 
probe it is mostly sensitive to the dipolar field from its 
nearest neighbour magnetic moment. Therefore, the x 
independent A, together with our previous conclusion re- 
garding the x independent size of Fe magnetic moment, 
imply that the implanted muons stop at a similar aver- 
age distance from a neighbouring Fe in all samples. How- 
ever, one expects that the implanted muons have higher 
probability to stop closer to a Fe ion for larger x. This 
contradiction may be explained if the muon stops within 
the van der Waals gap in the material, as we concluded 
from the small relaxation rates measured in the mother 
compound. In such a case the distribution of dipolar 
fields experienced by the muons from the Fe ions, which 
occupy a site in the Bi layer, does not depend stron gly 
on x. The width of this distribution can be calculated 23 
by distributing x/2 Fe ions randomly on a monolayer of 
a triangular lattice (of Bi, lattice constant a), and sum- 
ming up their contribution on muons which are located 
4.76 A away from it (within the van der Waals gap). Al- 
though the proposed model is very simple, we find that 
for an average Fe moment of 0.34/xb we obtain values of 
A which are consistent with the experimental results. 

The measurements presented here indicate that the 
low and high temperature magnetic phases differ in the 
underlying interactions. The high temperature regime 
seems to be associated with Fe-Fe coupling, and there- 
fore depends on x which affects the average distance be- 



tween Fe ions. In contrast, the weak x dependence at low 
temperature is strong evidence that the magnetization is 
dominated by the Fe single ion properties and its interac- 
tion with the Bi2Se3 lattice, e.g. crystal field etc. Hence, 
as the temperature is decreased the smaller single ion 
interaction with the lattice becomes more relevant, al- 
tering the high temperature size and ordering of the Fe 
moments. This is consistent with the strong anisotropy 
detected in magnetization measurement, which show a 
strong tendency of the Fe moments to align perpendicular 
to the c-axis, especially at low temperatures. Moreover, 
considering the low density of Fe ions, and therefore the 
fairly long range (a few a) magnetic Fe-Fe interaction, 
it is most probably mediated by charge carriers which 
extend over these length scales. Such effects have been 
observed in Fe doped Bi2Te3, though with magnetic tran- 
sition at much lower temperatures (< 15 K) 24 which also 
depend strongly on x. This was attributed to donor type 
Fe doping in both Bi2Te3 and Bi2Se3. However, in that 
study it was found that Bi 2 _ x Fe a; Se3 with x — 0.04 is 
non-magnetic. Furthermore, that the magnetic transi- 
tion temperatures in Bi 2 _ a; Fe a ;Te3 depends strongly on x 
and that the easy axis is along c. Finally, it is important 
to point out here that recently Lasia et alP3 predicted 
that in a magnetically doped TI the magnetic moments 
order ferromagnetically at zero tempertature, but udergo 
a first order transition to a spin density wave (SDW) at 
higher temperature before they reach the paramagnetic 
phase. Our results may be consistent with this model, 
though it is unlikely that the high temperature phase, 
which can be detected clearly in dc magnetization mea- 
surements, is due to a SDW ordering. Such ordering is 
typically associated with small bulk magnetization. 



magnetic at a relatively high temperature of ~ 250K. 
We identify two magnetic regimes; the first between 
~ 10 — 250K and the second below ~ 10K. Bulk magne- 
tization measurements show that the average size of Fe 
magnetic moment depends on temperatures but is almost 
independent of doping, and that the at the Fe moments 
align perpendicular to the crystallographic c-axis. Our 
/iSR results indicate that the low temperature regime is 
independent of doping level and is more likely determined 
by Fe single ion properties and crystal field effects. In 
contrast, the nature of the magnetism in the high temper- 
ature regime differs, such that Fe-Fe interactions become 
more important, and therefore the local magnetization 
depends more on doping. The observed magnetism is ho- 
mogeneous in the full volume of the sample and cannot 
be attributed to magnetic impurities. Finally, we con- 
clude that the appearance of local static magnetic field 
at high temperatures and low doping is strong evidence 
that the Fe-Fe interactions are dominated by conduction 
electrons, e.g. RKKY type interaction. Therefore, near 
the surface, we expect that the magnetization will be 
affected differently by the topological surface states in 
the two temperature regimes. At low temperatures, the 
broken crystal symmetry will affect crystal field, while 
the change in the nature of conduction carriers due to 
topological surface states will affect the high tempera- 
ture magnetism. Low energy /iSR measurements are un- 
derway to explore the magnetism and its dependence on 
surface proximity. The results presented here will pro- 
vide an excellent point of reference to detect any changes 
in the nature of the magnetism near the surface of Fe 
doped Bi 2 Se 3 . 



V. SUMMARY AND CONCLUSIONS 

In conclusion, we find that even at Fe doping levels 
as low as 5%, the full volume of the Bi2Se3 becomes 
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